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It was recently shown that dielectronic recombination measurements can be used for accurately 
inferring changes in the nuclear mean-square charge radii of highly-charged lithium-like neodymium 
[Brandau et al., Phys. Rev. Lett. 100 073201 (2008)]. To make use of this method to derive 
information about the nuclear charge distribution for other elements and isotopes, accurate electronic 
isotope shift parameters are required. In this work, we calculate and discuss the relativistic mass- 
and field-shift factors for the two 2s 2 5*1/2 — 2p 2 P°/ 2 3/2 transitions along the lithium isoelectronic 
sequence. Based on the multiconfiguration Dirac-Hartree-Fock method, the electron correlation and 
the Breit interaction are taken into account systematically. The analysis of the isotope shifts for 
these two transitions along the isoelectronic sequence demonstrates the importance and competition 
between the mass shifts and the field shifts. 
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I. INTRODUCTION 

Isotope shift (IS) measurements have a long tradition and are known as a valuable source for information about 
(changes in) the nuclear charge radii and distributions. Adding one or several neutrons makes the nuclei not only 
heavier, but may also modify its charge distribution quite remarkably During the last decades, optical laser 

spectroscopy considerably contributed to the field as an important alternative to a number of earlier x-ray measure- 
ments in muonic atoms or electron scattering experiments [H, Q . Using optical laser spectroscopy, measurements 
became possible for quite long sequences of stable and radioactive isotopes [3, H| . A particular merit of this method is 
that it provides nuclear-model-independent data if the electronic part of the isotope shifts is known sufficiently well. 
However, the determination of the electronic parameters through ab initio calculations is still a challenge for atomic 
theory owing to the high sensitivity of the IS parameters with regard to electron correlation and the representation of 
the wave functions, especially for (nearly-) neutral atoms and ions [9T4ll|. The electronic contributions to the isotope 
shifts can be theoretically attained with relatively high precision and reliability for few-electron ions, but as pointed 
out in Ref. fl2j , it would be still desirable to have more accurate results for the electronic parameters with up-to-date 
theoretical methods. 

Following the theoretical investigation on the effects of the nuclear charge distribution upon the total cross-section 
of dielectronic recombination for Li-like heavy ions by §chiopu et al. [l]|, isotope shifts of the 2s 2 £1/2 — 2p 2 Py 2 3 , 2 
transitions have been measured recently at the experimental storage ring at GSI Darmstadt (Germany) for the two 
stable even-even isotopes 142 Nd and 150 Nd of Li- like neodymium PJj]. From these measurements, it was shown that 
dielectronic recombination supplies a new technique to derive information on the nuclear charge distributions of 
heavy stable and unstable elements, provided the electronic parameters are available from theory. As far as Li-like 
ions are concerned, many investigations were performed. However, most of them are focused on neutral Li fl5j (and 
references therein) and Be + ion [2|, [HI, E3 where nuclear halo structures were predicted some years ago. Godefroid 
et al. have studied the isotope shifts of low-lying levels for Li I through O VI using the multiconfiguration Hartrec- 
Fock method [Hj]. Kozhedub et al. [l9| calculated recently the relativistic mass shift for the 2s 2 S , 1 / 2 — ^p 2 Py 2 an< ^ 
2s 2 S\/2 — 2p 2 P^ 2 transitions along the Li-like isoelectronic sequence. Nevertheless, the isotope shift parameters are 
scarce for Li-like ions, especially for the field shift electronic factors. 

In this work, we present the electronic parameters relevant to the isotope shift for the two 2s 2 S\/2 — %P 2 P°/2 3/2 
transitions in Li-like ions in the range from Z = 4 to Z = 90. The electron correlation and Breit interaction are 
treated within the framework of the multiconfiguration Dirac-Hartrec-Fock (MCDF) method. Comparison is made 
with other theoretical predictions, when available. Using our calculated parameters the balance between the mass 
shifts and the field shifts of these two transitions is discussed along the isoelectronic sequence with the assistance of the 
simple empirical formulae for the nuclear mass number and the nuclear mean-square radius. In addition we investigate 



in great detail the case of 150 - 142 ]Nfd 57+ for which experimental IS measurements were recently reported 14 



II. THEORY AND COMPUTATIONAL STRATEGY 
A. Isotope Shift Theory 

In the approximation of the first-order perturbation theory, the isotope shift, for a transition involving the atomic 
states u and I between two isotopes with mass numbers A > A' , can be parameterized as [20| 

8v A ' A ' = AK RMS M '~ M + F5(r 2 ) A > A ' , (1) 
MM' N ' w 

where AK RMS = (K RMS — K RMS ) refers to the relativistic mass shift parameter and F is the field shift factor; 
M' and M are the nuclear masses and S(r 2 ) A ' A represents the nuclear mean-square radius difference between two 
isotopes A and A'. 

The mass-shift parameter K is related to the expectation value of the relativistic recoil Hamiltonian 

K RMS = f(mRMs\*} ■ (2) 



The recoil Hamiltonian, correct to the order of (m/M) in the mass ratio of the electron and nucleus, is given by [2lM2 
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and is split into the one-body and the two-body relativistic mass shift operators 
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referred respectively as the relativistic normal mass shift (RNMS) and specific mass shift (RSMS) terms. Adopting 
this partition, the transition isotope shift ([1]) becomes 



8v A > A ' = (AK RNMS + AK RSMS ) — — + FS(r 2 ) A ' A ' 
V ; MM' V 1 



(6) 



Note that the normal mass shift values reported in the present work are estimated using the perturbation approach 
from the expectation values of the RNMS operator rather than adopting the scaling law (AK NMS = m e v) that 
becomes inappropriate for heavy relativistic systems [2ll . |25| . [26| . 

The field shift factor appearing in Eq. ([!]) has the form [27j 



3h \47reo 



A|*(0)| 



(7) 



which is proportional to the change of the electronic total probability density p e (0) at the origin between the two 
atomic states involved in the transition, i.e. A | ^(0)1 2 = p^(0) — pf{0), where p e (0) is defined by 



P e (0) = (*|^ ( 5(n)|*). 



(8) 



B. MCDF Method 

The multiconfiguration Dirac-Hartree-Fock (MCDF) method [28| is employed to yield atomic state wave functions 
(ASFs). On the basis of the MCDF method, the ASF is represented by a linear combination of configuration state 
functions (CSFs) with same parity P, total angular momentum J and its component along z-dircction Mj as 

^(PJMj) = c^iPJMj), (9) 

where {c^} are the mixing coefficients and {"/i} the sets of other quantum numbers needed for specifying the iV c 
CSFs. The latter are built from single-electron orbital wave functions. Applying the variational principle, the mixing 
coefficients Cj and single-electron orbital wave functions are optimized simultaneously via the self-consistent field (SCF) 
method. A more limited version is the relativistic configuration interaction (RCI) approach allowing only the mixing 
coefficients to be varied. In the RCI computation, the Breit interaction can be taken into account perturbatively as 
well. 

It should be emphasized that instead of a point-like nucleus, a finite nuclear charge distribution model, for instance, 
a uniform spherical model or a two-parameters Fermi model, must be used to generate the nuclear potential in the 
calculation in order to ensure the first-order perturbation to be valid for the field shift [29|. In the present work, the 
two-parameters Fermi nuclear model 

p(r) = ^-^T- (10) 

rw 1 _|_ e {r-c)/a ^ I 

is adopted. Here po is a normalization coefficient, c is the half-density radius and a = 41 * i3 is related to the surface 
thickness t of the charge distribution. In practice, the t = 2.30 fm value is used and c is computed according to the 



formulae given in Ref. [30( . Since the ASFs are insensitive with respect to details of the nuclear model [29j, |31( , the 
nuclear parameters from any stable isotope can be chosen to perform the SCF and RCI calculations. 
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TABLE I. Relativistic normal mass shift AK RNMS (in GHz u), specific mass shift AK RSMS (in GHz u) and field shift F (in 
MHz/fm 2 ) factors for the 2s 2 5*1/2 — 2p 2 P°/ 2 and 2s 2 Si/2 — 2p 2 Py 2 transitions in Be + . The label (n — x) specifies the 
orbital active set. 



2s 2 S 1/2 - 2p 2 P° j2 2s 2 S 1/2 - 2p 2 P° t 



Code 


Model 


^j^RNMS 




F 


AK RNMS 


^j>RSMS 


F 




DF 


-529.1 


-957.9 


-16.72 


-529.1 


-958.0 


-16.72 




SD(n=2) 


-521 


-2610 


-15.39 


-521 


-2610 


-15.39 


GRASP2K 


SD(n=7) 


-529 


-1028 


-16.94 


-529 


-1028 


-16.94 




SD(n=8) 


-538 


-1025 


-16.98 


-538 


-1025 


-16.98 




U SDT(n=5) 


-538 


-1025 


-16.99 


-538 


-1025 


-16.99 


MCDF-GME SDT(n=8) 


-532 


-1034 


-17.05 


-534 


-1034 


-17.15 



C. Computational Models 

The active space approach is utilized for building the ASF in order to capture the main correlation effects system- 
atically. In this framework, the configuration space is expanded by means of single (S) and double (D) excitations 
from occupied orbitals to some active set. The latter is augmented layer by layer, which makes it possible to monitor 
the convergence of the physical quantity concerned. Each correlation layer, labelled by n = 2, 3, 4, .. ., contains the 
s, p, d, ... orbitals (/max = 6). Starting with Dirac-Hartree-Fock (DF) calculations, MCDF calculations are carried 
out sequentially to optimize the correlation orbitals in the added layer, together with the mixing coefficients. The Brcit 
interaction is taken into account in the subsequent RCI computations. The calculations are performed by using the 
GRASP2K package [H, HH in which relativistic mass shifts corrections are implemented (24l |34| . For cross-checking 
our investigations on some cases, we also use the MCDF-GME package |35j . 



III. RESULTS AND DISCUSSION 

A. Electron correlation effects on isotope shifts 

Being interested in the entire isoelectronic lithium sequence and realizing the importance of correlation effects in 
the neutral end, we test our computational strategy on Be + . For this ion, we present in table U the relativistic 
normal mass shift (AK RNMS ), the relativistic specific mass shift (AK RSMS ) and the field shift (F) factors for the 
2s 2 Si/2 — 2p P°/ 2 an d 2s 2 5i/2 — 2p 2 Py 2 transitions, omitting the Breit interaction in the calculations. As can 
be seen from this table, SD excitations up to n = 8 are essential to bring these values to the demanded precision. 
However, one observes that the relativistic normal mass shift in these transitions is more sensitive to correlation due 
to large cancellation when making the difference between the level shifts considered. Complete convergence of these 
physical quantities is considerably hard to achieve, but we control the uncertainties within sa 1 % for the total mass 
shift and the field shift. 

Furthermore, we perform additional RCI calculations by appending the triple (T) excitations from the occupied 
orbitals to the active set of n = 5, to the largest SD configuration space (n = 8). The corresponding values marked 
as SDT are displayed in the same table. It is found that the contribution of triple excitations is definitely fractional. 
We should keep in mind that some orbitals are kept frozen in the layer-by-layer calculations to avoid convergence 
problems. The corresponding loss of orbital relaxation can be estimated with the MCDF-GME package [35| that 
allows to optimize all orbitals variationally. In this calculation, the SDT excitations up to n = 8 are used to generate 
the configuration space. The MCDF-GME results are listed in the last row in table U The satisfactory consistency 
between the results obtained by these two independent codes is on line with the 1% estimated uncertainty for the 
calculated shifts and demonstrates the adequacy of the computation strategy described in Section III CI that is used 
for heavier ions. 
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B. Electronic Factors in Isotope Shifts 

In table [TT1 we present the total mass shift AK RMS and field shift F factors for the 2s 2 Si/ 2 — 2p 2 -P 1 ° /2 an< ^ 
2s 2 Si/2 — 2p 2 Py 2 m Be + , Zn 27+ , Nd 57+ and Hg 77+ . In addition, the individual values of relativistic normal mass 

shifts AK RNMS and specific mass shifts AK RSMS are reported. To illustrate the effects of electron correlation and 
Breit interaction on these physical quantities along the isoelectronic sequence, the calculations are carried out using 
different computational models, namely, DF, MCDF, and MCDF+Breit. The differences between the MCDF and DF 
results reflect the correlations between the electrons. Since the correlation effects become usually less important for 
high-Z ions, we have reduced gradually the size of active orbital set with increasing the nuclear charge. Based on the 
MCDF model, the Breit interaction is included into the values denoted by the "MCDF+Breit" label. It is observed 
that both the electron correlation and the Breit interaction contribute little to the total mass shift and field shift 
factors in these two transitions, although the effect is larger for the AK RNMS contribution. Furthermore, we find out 
that the two-body part of the mass shift operator is at least two orders of magnitude larger than the one-body part 
in the middle- and high-Z region. This is in opposition with Seltzer j36| who pointed out that, in the case of K x-ray 
transitions, the specific mass shift is approximately —1/3 of the normal mass shift for 40 < Z < 70. 

We also compare our results with other available theoretical calculations [2 EMI that are collected in table HU As 
can be seen from this table, the agreement is satisfactory for the total mass shift AK RMS and the field shift F factors. 
However, we notice a quite large discrepancy with Kozhcdub ct al.'s value [l9[ in the RNMS of the 2s 2 S\/ 2 — 2p 2 P°i 2 
transition in Nd 57+ . Additionally, the differences between our Nd 57+ field-shift factors and Brandau ct al. results [14[ 
are relatively small (about 5%). Independent calculations [? ], carried out using the MCDF-GME package, are 
consistent with GRASP2K. 

The total results, obtained by the three computational models, for the relativistic mass shift and the field shift 
factors of the 2s 2 Si/ 2 — 2p 2 P 1 °/ 2 and 2s 2 Si/ 2 — 2p 2 Py 2 transitions are presented in table UTT1 and IIVI from Z = 4 

to Z = 90. It is worth noting that the trend of AK RMS for the 2s 2 Si/ 2 — 2p 2 P°p transition is not monotonous 
due to the fact that the relativistic corrections to the mass shift become important in the heavy elements. For 
this line, the contribution of the second term of Eq. © to the total line mass shift, counteracting the first term, 
increases drastically with the nuclear charge (10% in Fe 23+ , 47% in Nd 57+ and 90% in Th 87+ ). Note that the QED 
contribution to the recoil operator, calculated by Kozhedub et al. [l|| and omitted in the present work, becomes 
larger than the relativistic contribution for the high-Z region. The field shift factors for these two transitions however 
increase considerably rapidly along the isoelectronic sequence. 

The uncertainties of the mass- and field-shift parameters presented in table IIIII and IIVI mainly result from the 
electron correlation effects for the \ow-Z ions. As discussed in section 1111 Al the effects lead to about 1% errors along 
the isoelectronic sequence. Moreover, we neglect the frequency-dependent Breit interaction (FDBI) and the quantum 
clcctrodynamical (QED) corrections in the present calculations and these increase for the high-Z ions. We roughly 
estimate their contributions to the factors for Li-like Nd and Th ions. It is found that the FDBI and QED corrections 
affect the field shift factors by 1% for Li-like Nd and 3% for Li-like Th in both transitions under consideration. Also, 
they give 1% and 5% contributions to the mass shift factor in the 2s 2 S , 1 / 2 — 2p 2 P1/ 2 transition and 1% and 3% in 
the 2s 2 <Si/2 — 2p 2 Py 2 transition for Li-like Nd and Th ions, respectively. 



C. Mass and Field Shifts Balance 



The values of the mass shift and the field shift depend not only on the specific isotope pair, but on the transition 
under investigation as well. Consequently, the total isotope shift relies on the relative sign of the mass and field shifts 
if they are of the similar order of magnitude. In fact, the mass shift and the field shift are subject to different trends 
along the isoelectronic sequence, which makes it important to discuss the competition between these two expectation 
values in different Z-regions. 

To evaluate the mass shifts and the field shifts for the 2s 2 Si/ 2 — 2p 2 P°/ 2 and 2s 2 S , 1 / 2 — 2p 2 Py 2 transitions 
along the Li-like isoelectronic sequence, the difference of the nuclear mass and mean-square charge radius for a given 
isotope pair should be determined. For convenience, we use the mass number A instead of the nuclear mass and adopt 
the empirical formula [l3| 



1.98 + 0.015A 2 /3- v ' 

to inversely deduce A for a stable isotope of given atom with the atomic number Z. Considering two isotopes of the 
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TABLE II. Relativistic mass shift AK RMS (in GHz u) and field shift F (in MHz/fm 2 ) factors for the 2s 2 S 1/2 - 2p 2 P° /2 and 
2s 2 Si/2 — 2p 2 P%/2 in Be + , Zn 27+ , Nd o7+ and Hg 77+ . For comparison, individual relativistic normal mass shift /\j{ RNMS and 
specific mass shift coefficient AK RSMS (in GHz u) are presented as well. The label of n in parentheses represents the largest 
size of the active set. The number in square brackets is the power of 10. 



Model 




2s 2 Si/2 


- 2p 2 P° /2 






2s 2 Si/ 2 


" 2 P 2 P S/2 








Aif HMS 


F 






AK RMS 


F 










Be+ 










DF 


-529.1 


-957.9 


-1487 


-16.72 


-529.1 


-958.0 


-1487 


-16.72 


MCDF(n=8) 


-537.8 


-1025 


-1562 


-16.98 


-538.0 


-1025 


-1563 


-16.98 


MCDF+Breit 


-537.8 


-1025 


-1562 


-16.97 


-537.9 


-1025 


-1563 


-16.97 


Godefroid et al. [181 




-1035 




-16.91 




-1035 




-16.91 


Yan et al. [161 








-16.912 








-16.912 



DF -7475 -2177[2] 

MCDF(n=6) -7727 -2168[2] 

MCDF+Brcit -7895 -2167[2] 

Kozhedub et al. [19] -8054.2 -2165.449[2] 



Zn 27 + 

-2251 [2] -2009[2] 
-2245[2] -2012[2] 
-2246[2] -2009[2] 
-2246.00(3) [2] 



-1143[1] 
-1128[1] 
-1129[1] 
-1124. 97[1] 



-2194[2] 
-2187[2] 
-2188[2] 
-2168.232[2l 



-2308[2] 
-2300[2] 
-2301[2] 
-2300.73(3) [2] 



-2025[2] 
-2024[2] 
-2021[2] 



DF -1083[1] -8227[2] 

MCDF(n=5) -1342[1] -8196[2] 

MCDF+Breit -1449[1] -8196[2] 

Kozhedub et al. [19] -1641. 8[1] -8180. 90[2] 
Brandau et al. ^4| 
MCDF-GME 



Nd° 

-8336[2] 
-8331[2] 
-8341[2] 
-8345.08(25)[2] 



-7903 [3] 
-7929 [3] 
-7885[3] 

-7520[3] 
-7897[3] 



-8721[1] 
-8589[1] 
-8577[1] 
-8573. 3[1] 



-8768[2] 
-8761[2] 
-8775[2] 
-8769.29[2l 



-9640[2] 
-9620[2] 
-9632[2] 
-9626.62(25) [2] 



-8215[3] 
-8203[3] 
-8157[3] 

-7810[3] 
-8168[3] 



DF 

MCDF(n=4) 
MCDF+Breit 
Kozhedub et al. [19J 



4691[1] -1146[3] 
3990[1] -1140[3] 
2961[1] -1119[3] 



Hg"+ 

-1099[3] -5761[4] 
-1100[3] -5791[4] 
-1090[3] -5783[4] 
-1105. 6[3] 



-2193[2] -1452[3] 
-2158[2] -1451[3] 
-2181[2] -1434[3] 



-1671[3] -6239[4] 

-1667[3] -6225[4] 

-1652[3] -6213[4] 
-1669. 5[3] 



selected element, the S(r 2 ) A ' A can be obtained via the empirical formula (in fm) [37| 

(r 2 ) 1/2 = 0.836^1 1/3 + 0.570 (A > 9) . (12) 

Since the relativistic mass shift AK RMS and the field shift F factors are also Z-dependent, the mass shifts v MS 
and the field shifts v FS can be semi-quantitatively obtained as functions of the atomic number Z with the assistance 
of Eq. (fTTj) and (|T2^) . Using the MCDF+Breit results from table Hill and IIV1 we illustrate in figure [1] the absolute value 
of \v MS \ and \v FS \ for the 2s 2 S 1/2 ~ 2p 2 P° /2 and 2s 2 S 1/2 - 2p 2 P° /2 transitions in Li-like ions with 6 < Z < 90 
for the isotope pair (A, A + 1). The balance between the mass shifts and the field shifts can be easily observed from 
these two figures. As can be seen, in the region of 6 < Z < 35 the mass shift is somewhat larger than the field shift 
for the transitions concerned, but they are of the same order of magnitude. Therefore, both of them must be taken 
into account for a relevant analysis of isotope shifts. Over the isoelectronic sequence, the mass shifts can be regarded 
as constant, but the field shifts increase drastically towards high-Z values and become dominant at the end of the 
sequence. For instance, the field shift is one order of magnitude larger than the mass shift for Li-like Xc ion. As a 
result, one can safely neglect the mass shift for Z > 80 ions. 

In addition, it should be emphasized that the mass shift and the field shift cancel out with each other when building 
the isotope shifts of the 2s 2 Si/ 2 — 2p 2 P°j 2 and 2s 2 Si/ 2 — 2p 2 P^/ 2 transitions. Therefore, the total isotope shifts 
approach zero around Z = 34 (Se 31+ ). 



D. Isotope Shifts in iou ' 14/ NcT' + 

According to Eq. ([T]), we calculate the mass shift, the field shift and the total isotope shifts for the 2s 2 Si/ 2 —2p 2 P°/ 2 

and 2s 2 S\/ 2 — 2p 2 Py 2 transitions in the case of 150 > 142 Nd 57+ , using the MCDF+Breit mass- and field-shift factors. In 

this calculation, the atomic masses [38| are used to compute the mass difference between the nucleus 150 Nd and 142 Nd. 
The nuclear root-mean-square (rms) value (r 2 ) 1 / 2 taken from Ref. [l2j are adopted for estimating the corresponding 
<5(r 2 } 150,142 = 1.2909 fm 2 . Since Eq. (|T|) is based on the first-order perturbation theory, the corresponding results are 
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TABLE III. Mass- and field-shift factors, AK RMS (in GHz u) and F (in MHz/fm 2 ), for the 2s 2 S 1/2 - 2p 2 P° /2 transition for 
selected elements along the lithium isoelectronic sequence. The number in the square brackets represents the power of 10. 



AK RMS F 





DF 


MCDF 


MCDF+Breit 


DF 


MCDF 


MCDF+Breit 


"D„+ 


1 A O T roi 

— 1.487[o| 


— 1.5o2|3J 


— 1.562[3] 


1 Til 

— 1.d72[1J 


i ac\o \ 1 1 

— 1.698[1J 


— 1.o97[1J 




— 5. 540 [dj 


— 5.581[3] 


—5.581(3] 


— 1.408|2J 


1 A 1 o Tol 

— 1.418[2] 


— 1.417(2] 




i i nn r 1 

— 1.190[4] 


— 1.188 [4] 


— 1.188 [4] 


r ADA Tol 

— 5.434|2J 


— 5.454[2J 


—5.451(2] 


JNe 


— 2.050[4] 


— z.U4z|4J 


rt r\ a rt r a 1 

— z.U42[4] 


i ao a Tol 

— 1.486 |3J 


i /ion Tol 

— 1.489 [3J 


1 A o o Tol 

— 1.488 |3J 


c;ll+ 

Ol 


— 4.438[4J 


— 4.417 [4J 


— 4.41o[4] 


— 6.566 |3J 


—6.575 [3] 


— 6.571|3J 


A -15+ 

Ar T 


— 7.701 [4] 


— 7.668 [4J 


rr £?£?n r A 1 

— 7.669[4] 


i no o r a 1 

— 1.983 |4J 


i noc r a 1 

— 1.986 [4] 


i no A r A 1 

— 1.984 |4J 


Ti 19+ 


-1.182(5] 


-1.178(5] 


-1.178(5] 


-4.844(4] 


-4.850(4] 


-4.844(4] 


Fe 23+ 


-1.677 [5] 


-1.672 [5] 


-1.672(5] 


-1.033(5] 


-1.035(5] 


-1.033(5] 


Zn 27 + 


-2.251(5] 


-2.245 [5] 


-2.246(5] 


-2.009(5] 


-2.012(5] 


-2.009(5] 


Kr 33+ 


-3.252(5] 


-3.244(5] 


-3.245(5] 


-4.824(5] 


-4.833(5] 


-4.822(5] 


Mo 39 + 


-4.399(5] 


-4.392(5] 


-4.393(5] 


-1.052(6] 


-1.054(6] 


-1.051(6] 


Xe 51+ 


-6.997(5] 


-6.991 [5] 


-6.988(5] 


-4.176(6] 


-4.188(6] 


-4.167(6] 


Nd 57 + 


-8.336(5] 


-8.331(5] 


-8.341(5] 


-7.903(6] 


-7.929(6] 


-7.885(6] 


Yb 67 + 


-1.027(6] 


-1.027(6] 


-1.028(6] 


-2.165(7] 


-2.174(7] 


-2.172(7] 


Hg 77+ 


-1.099 [6] 


-1.100(6] 


-1.090(6] 


-5.761(7] 


-5.791(7] 


-5.783(7] 


Bi 80+ 


-1.072 [6] 


-1.073 [6] 


-1.053(6] 


-7.711(7] 


-7.752(7] 


-7.741(7] 


Fr 84+ 


-9.772 [5] 


-9.779(5] 


-9.379(5] 


-1.133(8] 


-1.139(8] 


-1.137(8] 


Xh 87+ 


-8.465(5] 


-8.472(5] 


-7.832(5] 


-1.512(8] 


-1.521(8] 


-1.518(8] 



TABLE IV. Mass- and field-shift factors, AK RMS (in GHz u) and F (in MHz/fm 2 ), for the 2s 2 S 1/2 - 2p 2 P 3 ° /2 transition for 
selected elements along the lithium isoelectronic sequence. The number in the square brackets represents the power of 10. 



AK RMS F 





DF 


MCDF 


MCDF+Breit 


DF 


MCDF 


MCDF+Breit 


Be+ 


-1.487(3] 


-1.563(3] 


-1.563(3] 


-1.672(1] 


-1.698(1] 


-1.697(1] 


C 3+ 


-5.543(3] 


-5.583(3] 


-5.584(3] 


-1.409(2] 


-1.418(2] 


-1.417(2] 


Q 5+ 


-1.191(4] 


-1.189(4] 


-1.189(4] 


-5.436(2] 


-5.453(2] 


-5.451(2] 


Ne 7+ 


-2.054(4] 


-2.046(4] 


-2.046(4] 


-1.487(3] 


-1.489(3] 


-1.489(3] 


Si n + 


-4.457(4] 


-4.435(4] 


-4.434(4] 


-6.576(3] 


-6.579(3] 


-6.575(3] 


Ar 15 + 


-7.761(4] 


-7.726(4] 


-7.727(4] 


-1.989(4] 


-1.989(4] 


-1.987(4] 


rpjl9 + 


-1.197(5] 


-1.192(5] 


-1.192(5] 


-4.864(4] 


-4.863(4] 


-4.857(4] 


Fe 23+ 


-1.707(5] 


-1.701(5] 


-1.701(5] 


-1.039(5] 


-1.039(5] 


-1.037(5] 


Zn 27 + 


-2.308(5] 


-2.300(5] 


-2.301(5] 


-2.025(5] 


-2.024(5] 


-2.021(5] 


Kr 33+ 


-3.378(5] 


-3.368(5] 


-3.369(5] 


-4.883(5] 


-4.880(5] 


-4.869(5] 


Mo 39+ 


-4.649(5] 


-4.638(5] 


-4.640(5] 


-1.070(6] 


-1.069(6] 


-1.065(6] 


Xe 51+ 


-7.787(5] 


-7.771(5] 


-7.771(5] 


-4.304(6] 


-4.298(6] 


-4.278(6] 


Nd 57 + 


-9.640(5] 


-9.620(5] 


-9.632(5] 


-8.215(6] 


-8.203(6] 


-8.157(6] 


Yb 67 + 


-1.307(6] 


-1.305(6] 


-1.306(6] 


-2.291(7] 


-2.286(7] 


-2.283(7] 


Hg 77+ 


-1.671(6] 


-1.667(6] 


-1.652(6] 


-6.239(7] 


-6.225(7] 


-6.213(7] 


Bi 80+ 


-1.775(6] 


-1.770(6] 


-1.745(6] 


-8.422(7] 


-8.401(7] 


-8.383(7] 


Fr 84+ 


-1.900(6] 


-1.894(6] 


-1.846(6] 


-1.253(8] 


-1.250(8] 


-1.247(8] 


Th 87+ 


-1.977(6] 


-1.970(6] 


-1.896(6] 


-1.690(8] 


-1.685(8] 


-1.681(8] 
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FIG. 1. Absolute value of the mass shifts \5v | (in blue) and the field shifts \5v F | (in red) as functions of the nuclear charge 
6 < Z < 90. All values in GHz. Results are shown for 2s 2 S\/2 — 2p 2 -P°/ 2 (upper figure) and 2s 2 S\/2 — 2p 2 -P3/2 (lower 
figure), calculated with correlated wave functions (MCDF+Breit) for the isotope pair (A, A+ 1). 



labelled as "PT" . The total isotope shifts for these two transitions are obtained by adding the mass shift and the PT 
field shift contributions. The experimental measurements [l4j for the isotope shifts are also displayed in this table. 
Considering that the contributions from the nuclear deformatio n lipp, n uclear polarization fl3 | , nuclear size effects on 
the one-loop QED corrections [3!| and QED recoil corrections |l4l. Il9l. are not taken into account in the present 
work, they are subtracted from the experimental values to comply with an objective comparison. The corresponding 
values are reported in the last line of Table |V] and are marked by a "*" . It is found from this table that the mass 
shift values are in perfect agreement with Brandau et al. [l4| and Kozhedub et al. results, and the isotope shifts 
are in satisfactory consistence with observations flil ]. 

One should keep in mind however that the 5(r 2 ) 150 ' 142 value of Angeli has been recently revised by Brandau et 
al. from a combination analysis Adopting this latest value, ie. S(r 2 ) 150 ' 142 = 1.36(1)(3) fm 2 for estimating the 
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TABLE V. Comparison of calculated isotope shifts and parameters for the 2s 2 Si/2 — 2p 2 P°/2 an d 2s 2 Si 



/a 



2p 



2 no 
^3/2 



transitions of lithium-like 150 ' 142 Nd j7+ with experiment and other computations. All shifts are given in meV. The nuclear rms 
radii of these two isotopes are taken from Ref. 





2s a S 1/a - 2p ~ z P° /2 


2s a S 1/a 


- 2p ^P 3 ° /2 




This work Others 114] 


This work 


Others 114] 


Mass shift 

S(Hrms) 


1.30 1.30 


1.50 


1.50 


Field shift 

PT 

VA 


-42.09 
-40.30 


-43.55 
-41.70 





Isotope shift 

PT 

VA 



Exp. ]14] 



-40.79 
-39.00 
-41.18 f 



-42.05 
-40.20 
-42.45 1 " 



-40.2(3)(6) 
-41.33* 



-42.3(12)(20) 
-43.4* 



2\ 150,142 



from 



instead of [1 



^VA values rescaled using the 5(r 2 ) 
* Obtained by subtracting the nuclear deformatio n 1391 ] . nuclear polarization |14|], QED recoil [14j and nuclear size effects on 
the one-loop QED corrections [3!| from the experimental values of Brandau et al. [Til ]. 



PT theoretical field shift leads to 5 % discrepancies in the isotope shifts. With respect to the accuracy of calculated 
field shift factors, it is worthwhile to point out that the higher-order nuclear moments (5{r 4 ), <5(r 6 ), . . . , etc.) [36[ are 
neglected in the PT field shifts. These missing contributions can be retrieved by estimating the FS from the difference 
of the transition energies calculated using the specific nuclear potential for each isotope. Adopting the rms values 
of Angeli [TJJ to build the nuclear potential of both isotopes 142 Nd and 150 Nd, we recalculate the field shifts in the 
MCDF+Breit model. The corresponding results are marked "variational" (VA) in table [V] A difference of ~ 5 % is 
found between the PT and VA field shifts, which accounts for the higher-order moments correction. Rcscaling our 
variational field shift with the latest mean-square charge radius difference value gives (5i? 150 ' 142 = —41.18 meV and 
—42.25 meV for the 2s 2 Si/ 2 ~ 2p 2 Py 2 an( l 2s 2 Si/ 2 — %P 2 Py 2 transitions, respectively, in very nice agreement 
with Brandau et al.'s measurements, when taking the nuclear deformation, QED recoil and nuclear contributions into 
account. 



IV. CONCLUSIONS 



Based on the MCDF method, we calculate the mass shift AK and the field shift F factors for the 2s 2 S 1/2 -2p 2 P° /2 
and 2s 2 S , 1 / 2 — 2p 2 Py 2 transitions in Li-like ions with 4 < Z < 90. The uncertainties of these parameters are 
controlled within 1% in the Breit approximation thanks to the rather complete description of electronic correlation. 
Our results provide the reader with a fast estimate of the expected IS for a given element and pair of isotopes, using 
formula (fT]). If the isotope shift is measured for one of the lines, alternatively, the AK RMS - and F-parameters can 
be utilized to extract the change in the mean-square charge radii. 

Using these results we analyze the competition between the mass shift and the field shift for the 2s — 2p resonance 
doublet along the isoelectronic sequence. It is found that the mass shift and the field shift possess similar orders of 
magnitude in the Z < 40 range, so that one should consider both of them for a relevant analysis of isotope shifts, 
especially for extracting the nuclear mean-square charge radius. Towards the high-Z region, the field shift grows 
rapidly and becomes dominant at the end of the sequence. Therefore, the mass shift can be safely neglected in the 
region of Z > 80. In addition, the total isotope shifts approach zero around Z = 35 due to strong cancellation between 
the mass and the field shift contributions. 

A detail analysis of the isotope shifts in the case of 150 ' 142 Nd 57+ shows that higher-order nuclear moments, often 
neglected in the calculation of the field shift, should be considered for very highly charged ions in order to extract 
5(r 2 ) values from experiments. 
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